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Abstract 

The conductance of systems containing two tunnel point-contacts and a single subsurface scat- 
terer is investigated theoretically. The problem is solved in the approximation of s-wave scattering 
giving analytical expressions for the wave functions and for the conductance of the system. Con- 
ductance oscillations resulting from the interference of electron waves passing through different 
contacts and their interference with the waves scattered by the defect are analyzed. The prospect 
for determining the depth of the impurity below the metal surface by using the dependence of 
the conductance as a function of the distance between the contacts is discussed. It is shown that 
the application of an external magnetic field results in Aharonov-Bohm type oscillations in the 
conductance, the period of which allows detection of the depth of the defect in a double tip STM 
experiment. 

PACS numbers: 61. 72. J- Point defects and defect clusters; 73.63.Rt Nanoscale contact; 74.55.+V Tunneling 
phenomena: single particle tunneling and STM 
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With the further development of scanning tunnelling microscopy (STM) it has become 
clear that a single STM-probe is often not enough for obtaining information on the detailed 
characteristics of the surface under investigation. A logical development of the one-tip 
approach is a dual-tip experimental setup, which can provide us with richer information 
than conventional single-probe STM. Despite the apparent technical complexity of the dual- 
tip STM (DSTM) in comparison with standard STM several groups have demonstrated 
successful solutions for such refinement of the STM-technology [lH4|. 

DSTM can be realized in different ways. For example, it can be a spatially extended STM 
tip with two protrusions, each ending in a cluster or a single atom [5]. A second approach 
is a coaxial beetle-type double-tip STM design that looks advantageous in retaining the 
standard STM stability [§J. The most versatile DSTM comprises two individual tips, which 
can be driven independently. In this case the distance between the tips is limited in principle 
only by a parameter such as the characteristic tip radius [2]. Another original example of 
the DSTM was proposed in where one contact can be created directly on the surface, 
while the other one was the STM-tip itself. 

For DSTM experiments with two independent probes there are different possibilities 
for applying voltages to the tunnelling contacts. There are two basic circuit designs: In 
the first one electrons are emitted from the first contact and then gathered at the second, 
i.e. the current flows from one contact to the other through the surface being probed 

0]. This method allows capturing a trans-conductance map, and in addition allows the 
implementation of three-terminal ballistic electron emission spectroscopy (BEES) without 
introduction of macroscopically bounded contacts [3J. In the second basic scheme proposed 
in Ref. [BJ the bias is applied between the two tips and the sample, i.e. the current flows 
from two contacts into the sample. 

Subsurface defects, adatoms, and steps on the metal surface result in the appearance of 
Friedel-like oscillations in the STM conductance G = dl/dV - a nonmonotonic dependence 
of G with the distance between the STM tip and the defect r$ (for a review see [10(). The 
study of this dependence can be used for the detection of buried defects and for investigation 
of their characteristics. Methods for determining defect positions below a metal surface 
using a single tip STM have been proposed before: this can be achieved using the period of 



oscillation of the conductance as a function of bias [11|, [12J or by exploiting the interference 



pattern of conductance as a function of position, G (r ), which is very pronounced for open 



directions of Fermi surface |13Hl5l]. These approaches are very suitable for the surfaces 
of simple metals, such as the noble metals, but application to conductors having a more 
complicated Fermi surface geometries will be difficult and has not yet been explored. 

In the present work we examine the case of injection of electrons to the surface by the 
first and the second contacts simultaneously. We consider this realization of a double-tip 
experiment as a natural refinement of the sin gle- tip STM problem for the study of single 



defects buried under the metal surface [ll|, 112, [15 . 

The idea of using multiple tunnelling contacts for determining the depth and location 
of impurities under a metal or semiconductor surface has been expressed earlier in Ref. jl] . 
The paper by Niu et al, Ref. proposes a method for determining the desired depth 
by measuring the trans-conductance between two tips of the dual-tip scanning tunneling 
microscope. In the present paper we propose a different approach, namely, by measuring 
the phase change Atf in the conductance oscillations as a function of the distance between 
two STM tips d. Such phase changes can be measured experimentally with great precision. 
We show that A$ can be expressed in terms of the distance d (in units of the Fermi wave 
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FIG. 1: Schematic arrangement of the system of two tunnel contacts, modelled as two orifices in an 
infinitely thin interface between conducting half-spaces. The inset shows the equivalent circuit with 
two STM tips, which provide the electron tunnelling paths through small areas with characteristic 
radii a% and 02- 



vector /cf), the position of the defect in the plane parallel to surface plane p , which is easily 
denned experimentally, and the unknown depth of the defect zq ■ Thus by measurement of 
Ai? (d) it is possible to determine the depth of the buried impurity. The procedure of defining 
Zq is further simplified when a magnetic field H is applied to the system. In this case the STM 
conductance G undergoes Aharonov-Bohm type oscillations. These oscillations result from 
the quantization of the magnetic flux through the area formed by the electron trajectories 
from the contacts to the defect and the line connecting the contacts (FigJT]). For a weak 
magnetic field the electron trajectories and the line connecting the contacts form a triangle, 
and from its area S the defect depth can be found easily. 

As a model for the double-tip STM geometry we consider two metal half-spaces separated 
by an infinitely thin nonconducting interface at z = 0, which contains two small regions 
(contacts) that allow electron tunnelling (see FigJT]). The origin of the coordinate system 
r = is chosen in the center of the first contact. The x axis is directed along the the line 
connecting the contacts. For the potential barrier in the plane z = we use the function 
17 

U(r) = U f(p)5(z). (1) 

In our case / (p) describes two "windows" for electron tunneling and the reciprocal function 
f^ 1 (p) can be presented as a sum of two terms 

/ _1 (p) = x(p/ a i) + x(|d-p|/a 2 ), 2tc J dxxx(x) = 1, (2) 



where % ( x ) — 1 f° r x ~ 1> an d X ( x ) ^ 1 f° r x ^> 1, a\ t 2 are the characteristic radii of 
the contacts, p is the component of the vector r parallel to the plane z = 0, d is a two 
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dimensional radius vector from the center of first contact to the center of second one. The 
absolute value d is the distance between contacts, assuming that this is smaller than the 
shortest relaxation length. 

In the vicinity of the contacts a single defect is placed described a short range potential 
D(t), 

D{T) = gD {\T-r \), (3) 

where g is the constant of interaction of the electrons with the defect, and D Q (\r — r |) is a 
spherically symmetric function localized within a region of characteristic radius centered 
at the point r = r , which satisfies the normalization condition 

oo 

4vr J dr'r' 2 D (r') = 1. (4) 
o 

For calculation of the conductance G we proceed as before. The probability density 
current is found by using the wave function ip (r) for the electrons tunnelling through the 
potential barrier in the plane of the orifices. The total electric current / in the system is 
calculated by integrations over electron momenta and over a real-space surface overlapping 
the contacts. We will take the temperature to be zero, and assume a small applied voltage 
V such that we stay in the linear regime of Ohm's law, / = GV . Under these assumptions 
the conductance G can be written as 

G = (e F ) f dn p [ dtlr 2 Im ty>* (r) (r)] . (5) 

TJX j J 

S F ,v z >0 S 

In Eq.flH]) m* is the effective electron mass, v (sp) is the electron density of states at the Fermi 
level, dQ and dfl p are solid angles in the real and momentum spaces, respectively. As the 
surface for space integration we choose a half-sphere of radius r, larger than distance between 
the contacts d and centered at the center of first contact, r = 0, and covering the contacts in 
the lower half-space, z > 0. The integration over the directions of the momentum over the 
Fermi surface Sp is carried out for electrons tunnelling and having a positive projection v z 
of the electron velocity on the contact axis z. As a consequence of the conservation of total 
current the integral over dQ does not depend on the length we choose for the radius r. 
The electron wave function ip (r) satisfies the Schrodinger equation 



V 2 + ^(,-D(r)) 



^W = 0, (6) 



subject to the boundary conditions of continuity and of the jump of its derivative at z = 0. 
In Ref.[l7| a solution of Eq.()6]) was found for an arbitrary function f (p), in the limit of 
weak tunnelling, 1/Uq —> 0, and for a purely ballistic contact (no defects present), 
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% (P, z) = -— ?477 / dVe^V^ / dp' 6 —— , (7) 
(2tt) m*U J J f (p'J 



where k' z = a/x 2 + k 2 — h' 2 , and x and k z are the components of the vector k parallel 
and perpendicular to the interface, respectively. As a special case the authors of Ref.[l7| 
considered a system of several orifices with different radii. 
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FIG. 2: Squared modulus of the wave function (|10p . The defect sits at kr® 
distance between the contacts kd = 14, the scattering phase shift is Sq = 1.5. 



(10,0,15), the 



The characteristic radius of the region through which the electrons tunnel from the STM 
tip into the sample has sub-atomic size (a ~ 0.L4) while the Fermi wave vector is fcp ~ 1A _1 . 
By using the condition kpai^ < 1 we find, after integrating over x' in Eq.(j7j), 



it(k z 



{ka{f -hfikr) + {ka 2 f -,h? (kr')e i>ed 



where h^{kr) is the spherical Bessel function of the first order, r' — |r — d| , and 



t(k z 



h 2 k z 

im* Uq ' 



(8) 



(9) 



is the transmission amplitude of the electron wave function passing through a homogeneous 
barrier. Note that in the limit kpai 2 — > the result of Eq. (JSJ) does not depend on the 
concrete form of the function x{p/ a ) m Eq. (j21) and the wave function Eq. (jSJ) as well as 
the conductance of the system are expressed in terms of the effective areas of the contacts, 



na 



1,2- 



The effect of electron scatte ring by the short-range potential can be taken into account 
by the method proposed in Ref. 18(. If the radius of action of the potential D(r) is of the 
order of Fermi wave length Ap, in the region of the defect |r — r | < rjj the wave function 
if) (r) can be taken as a constant if) (r ) , as for a 5-function. Under this approximation 
Eq. takes the form of a non-homogeneous equation with the right-hand member being 
-^-Dijj (r ) . In the limit 1/Uq — > a solution of this equation can be expressed in terms 
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of the solution of the homogeneous equation (see Eqs. (JTj), flB])) and the retarded electron 
Green's function of Eq.Q for the semi- infinite half-space 

^ (r) = % (r) + 2 ~^T (k) [G« (|r - r |) - G« (|r - r„|)] Vo (r„) , (10) 

where ro = (p , —zq). T (k) is the scattering matrix, which for a short-range scatterer can 
be expressed in terms of the s-wave scattering phase shift So 2Q 



T (A;) = 7 ^ v . (11) 

The Green's function 

G« (*) = - 6XP 4 ^ X) (12) 

is the retarded Green's function of a free electron. The phase shift 5q is determined by the 
scattering strength g as, 

(oo 
l-^f<' W D H ] . (13) 

Figf2] illustrates the spacial variation of the wave function ( flOj) for the case when the 
contacts and the scatterer are all placed in the plane y = 0. The interference of election 
waves passing through different contacts and their interference with the waves scattered by 
the defect are clearly visible. In order to make the effects more visible we used in Figj2] a 
large value for the scattering phase 5q = 1.5, which is acceptable only for Kondo resonance 
scattering by a magnetic impurity (see, for example, jl9j]). The grey circle round the point 
r = r in Figj2]is the region in which the Eq.f jT0|) is not valid because of divergence of the 
Green function. 

Substituting the wave function f fTOj) into the general expression for the conductance G 
flS} we find 

G = G c (ai, a 2 ) + G osc , (14) 
where G c is the conductance of the double contact system in the absence of the defect 

G c (oi, a 2 ) = G (oi) + G (a 2 ) + G 12 , (15) 

Go is an inherent conductance of the single contact [Pi 

Go(a) = |t(A; F )| 2£! |^, (16) 

and G12 takes into account the interference of electron waves passing through different 
contacts 

G 12 = \t{k F )\ ^-^ / (k F d). (17) 

Here we introduced the notation 

= (18) 
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ji (x) the spherical Bessel function of the first kind such that / (0) = 1. The second term 
in Eq. (ll4p . G osc , describes the quantum interference resulting from the scattering of the 
electrons by the defect 

G osc (r , d) = G K) r (r ) + G (a 2 ) T (r ) + 2G 12 ^ (r , r' ) // (k F d) , (19) 

where r = (p — d, z ), and r' is the distance between the defect and the second contact. 
The functions T (r ) and T (r' Q ) take into account the effect of interference of electron waves 
passing through the contact and returning to the same contact after scattering by the defect, 

r (r) = — L sin 5 ^ [Uj, (k F r) 7 (k F r) + 6(1- j (2k F z)) {k F r)~ A {{krf + l) sin 5 ] , 

(20) 

where 



F(z) = l + 2sin<5 

and 



2 (2k F z) 2 



- jo {2k F z) sin 5 - y (2k F z) cos S Q 



(21) 



7 (r) = -y x (k F r) cos 5 + sin S {j x (k F r) (j (2k F z) - 1) + y (2k F z) y x (k F r)} . (22) 

In the last term in Eq. (|T9l \l/ (ro, r ) describes the interference of electron waves that arrive 
at the other contact after scattering by the defect, 

* ( r , r') = F- 1 sin <5 ^ [ji {k F r') 7 (r) + j, (k F r) 7 (r') + (23) 

+ sin 5 (jo (k F d) - j (kW^ + d 2 ) j x 
(ji (^) ji (^') - V\ (k F r) y x {k F r'))\ , 

For a 2 = (i.e., when we have just a single contact) Eq.f ll4p coincides with the expression 
for the conductance of a tunnel point contact obtained in Ref.|20]. Fig. |3] illustrates the 
dependence of the oscillatory part of the conductance (fT9|) on the position of the defect in 
the plane z = z Q . The oscillatory pattern presented in Fig. |3] represents an image which 
could be obtained by DSTM when mapping the tunnelling conductance in the vicinity of 
the subsurface defect. 

The general formula for the conductance ( !T4l can be simplified for large distances be- 
tween the contacts and the defect, r ,r' Q ^> l/k F , and for a weak scattering potential 
5q ~ —gm*k F /27rh 2 1. Under these assumptions the normalized oscillatory part of 
the conductance, in the linear approximation in g, can be written as 

^ = -^o4 ( \ sin 2k F r + -1 sin 2k F r> - sin k F (r + r' )\, (24) 

For simplicity we take here a\ = a 2 = a. Equation (124"]) shows that in contrast to one tunnel 
point contact, for which G osc ~ sin2/c F ^o when k F r ^> 1, the oscillatory dependence of the 
double contact has a phase shift $ that depends on the distance between the contacts 

^£ „ sin (2A; F ro + i?) , (25) 
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FIG. 3: Dependence of the normalized oscillatory part of the conductance G osc /Gq as a function 
of the defect position p Q in the plane parallel to interface z = 0, kpzo = 5 . The distance between 
the contacts is taken as kpd = 20, and the scattering phase shift is So = 1.5. 



where 



$(p , z o, d) — — arcsin 



sin 2ip + 2f(kpd) sin (p 



'2 + Af 2 (k F d) + 4f(k F d)(cos 2<p + 2 cos y) 



d <C r , 



(26) 



and Lp = kF ^° d . The defect position in the plane parallel to the surface, p , is known from 
the interference pattern of the conductance oscillations (see FigOJ). In principle, the depth 
of the defect zq may be found from the experimental data in the following way: changing the 
distance between the contacts over a small range Ad <C d leads to the appearance of an ad- 
ditional phase shift A$ (p , z , d), which can be defined from the dependence G osc (p , z , d), 
see Fig. fflfy. The depth z can be obtained as a numerical solution of the equation 



A7? = 1 ? , d (p 0j z 0j d) Ad. 



(27) 



Let us now consider applying a magnetic field H parallel to the surface of the sample (see 
Fig. [1]). If the external magnetic field is sufficiently weak, such that the radius of the electron 
trajectories r# = hckp/eH is much larger than the distances between the the contacts and 
the impurity, ro, r' Q , the magnetic distortions of the trajectories j2l| are negligible, i.e. the 
trajectories can be considered as straight lines. 

Under this condition of r H 3> r , r' , the zero-field wave- function tp (r) acquires an addi- 
tional phase: 
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FIG. 4: Dependencies of the oscillatory part of the conductance on the coordinate x$ of the defect 
for different distances between the contacts, for yo = 0, k F zo = 5, and 5q = 0.1. 



te 



i/j{r) = V(r)exp | / A(i" ) dr" 



(28) 



and the Green function similarly takes the form [16 



le 



G (r, r ) = G (r, r ) exp ( - / A (r") dv" 



(29) 



ro 



Here, A (r) is the vector potential of the magnetic field. 

On account of this change in the wave function (1281) the formula for the conductance G 
is modified and takes the form: 



G = G c (ai, a 2 ) + G osc (r , d, H) 



(30) 



G osc (r , d, H) = G (a x ) T (r ) + G (a 2 ) T (r' ) + 2G 12 V (r , r' ) // (k F d) , (31) 



# (r, r'l 



1 • A Z 
— smo — : 



7T$ 7r$ 
ji (fepr') ( 7 (r) cos — + 7 (r) sin — 

^0 ^0 



+ 



(32) 



7T$ 71$ 

+Ji (for) ( 7 (r ) cos — + 7 (r ) sin — 
<±> 



+ 



+ sin 5 Q (j (k F d) - j (k F V4:Z 2 + d 2 j J x 
(jl (kpr)ji (k F r') - y 1 (k F r) y 1 (k F r'))} 
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7 (r) = cos5 ji (for) + smS {yi (for) (j (2k F z) - 1) - y (2foz) jt (for)} , (33) 

where 7 (r) is denned by ( 122]) . $ = ffch/e is the flux quantum and $ = HS is the magnetic 
flux through the triangle formed by vectors r , r , and the vector d connecting the contacts. 
At H = the expression (I3"U|) reduces to the formula obtained earlier (TH|) . 
At r , r ^> 1/&F, and 5 <C 1, Eq.f l3Tj) takes the form 

AG 05C (r ,rf,ff) ^ 12fr)zg/ (fc F rf) ^ 
G k F (r r' ) 2 

cos k F r Q sin /c F r n + cos k F r n sm k F r Q 

V $0/ \ $0/ 



Similar oscillations in the electron local density of states have been predicted in Ref. [22 
for a system of two adatoms and an STM tip in a plane perpendicular to a surface magnetic 
field. ' 

If the period of the oscillations is known, the depth z can be determined using the 
following procedure: In the most convenient geometry of the experiment the contacts should 
be placed so that the vectors ro, r and the normal to the sample surface are situated in the 
same plane, i.e. the vectors H and S are parallel. For our illustration in FigJTJthat means 
the coordinate p of the defect in the plane xy is on the line connecting the tips. In this 
case the relation between the period of oscillations AH and the depth zq is very simple 

Zo = dKH (35) 

Note that observation of the conductance oscillations AMI) requires a sufficiently strong 
magnetic field. Currently in low temperature STM the magnetic field up to 15T is reachable 



23j |. For example, in order to observe the quote of period AH for z = d =20 nm it is 
necessary to apply the field H = 5T. For typical metals, for which X F ~ 0.1 nm, for the 
distance between the contacts and the defect r > 10 nm the amplitude of conductance 
oscillations become very small G 

osc rsjl GUAp/ro) ~ (10~ 4 ^ 10- 5 )G - Therefore more 
suitable objects for application of proposed magnetic method of determination of the defect 
position below the surface are semiconductors, semimetals (Bi, Sb and their ordered alloys) 
where the Fermi wave length Af ~ 10 nm. Also the large amplitude G osc ~ (10~ 2 -f- 10 -3 ) Go 
could be expected in the metals of the first group, a Fermi surface of which has small cavities 
with effective mass m* ~ 10~ 2 -j- 10~ 3 mo (mo is the mass of a free electron). As well a 
low temperature STM should be used to avoid electron-phonon scattering on the electron 
trajectory. 

Thus, in this paper we have investigated theoretically the conductance of the system 
consisting of two close tunnel point contacts in the vicinity of which the point defect is 
situated. In approximation of s— wave scattering which is valid for short range scattering 
potential the oscillatory dependence of conductance on the separation between the contacts 
and their distances from the defect is studied. We proposed an alternative way that allows 
to determine the depth of the subsurface impurity by measuring the phase change in the 
conductance oscillation, arising when we change the distance between the contacts of double- 
tip STM. Also it was obtained that when the case of low magnetic field which is parallel to 
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the surface of the sample the depth of the subsurface impurity can be easily found from the 
period of Aharonov-Bohm type oscillations of conductance, which arise in this case. 
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